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Abstract 

We show from a simulations-based study of the TACTIC telescope that fractal and wavelet 
analysis of Cerenkov images, recorded in a single imaging Cerenkov telescope, enables al- 
most complete segregation of isotropic gamma-ray initiated events from the overwhelming 
background of cosmic -ray hadron-initiated events. This presents a new method for mea- 
suring galactic and extragalactic gamma-ray background above 1 TeV energy. Preliminary 
results based on this method are reported here.Primary aim is to explore the possibility of 
using data recorded by a single imaging atmospheric Cerenkov telescope(.L4CT) for mak- 
ing accurate measurements of diffuse galactic and extragalactic gamma-ray flux above ~ 
1 TeV energy.Using simulated data of atmospheric Cerenkov images recorded in an IACT, 
initiated , both, by cosmic ray protons and diffuse gamma-rays with energies above 4 TeV 
and 2 TeV respectively, we identify the most efficient fractal/wavelet parameters of the 
recorded images for primary identification. The method is based on the pattern recognition 
technique and employs fractal and wavelet analysis of the recorded Cerenkov images for 
gamma-hadron segregation. We show that the value of wavelet dimension parameter Bq can 
segregate Cerenkov images of hadronic origin from those of diffuse gamma-ray origin with 
almost 100% accuracy. We use the results to get a preliminary upper limit estimate of the 
diffuse galactic gamma-ray flux within galactic range of |6|< -5° and \l\< 200° above 2 
TeV from a 36h data set recorded by the TACTIC telescope. 
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1 Introduction 

The diffuse high energy gamma-ray background from the galactic plane, mapped 
extensively at energies above 100 MeV by the EGRET experiment, is generally 
attributed to nucleon-nucleon scattering and the decay of neutral pions produced 
in cosmic ray interactions with the interstellar matter (ISM) [1-5]. As the galactic 
cosmic ray spectrum, from sources like SNR and pulsars,is expected to extend upto 
~ 10 15 eV, the diffuse gamma-ray spectrum resulting from cosmic ray interactions 
with the ISM can easily extend into the 10 s of TeV range. Imaging atmospheric 
Cerenkov telescopes, which operate at energies of ~ 1 TeV, can look for galactic 
diffuse gamma-ray emission especially if two or more such systems are operated 
in the stereoscopic mode to enable accurate primary arrival direction reconstruc- 
tion. However, imaging Cerenkov telescopes have an inherently low sensitivity 
for differentiating between background cosmic-ray initiated events and the off- 
axis(diffuse) gamma-ray initiated events because the image shape and orientation 
parameters only separate primary gamma ray-initiated events with arrival direction 
along the principal axis of the light collector[6] . There have only been a few reports 
of searches for diffuse galactic gamma-ray emission using imaging Cerenkov tele- 
scopes and in several cases upper limits on the flux at energies > 1 TeV have been 
placed [7-10]. The first unambiguous detection of > ITeV diffuse gamma-rays 
from the galactic plane has recently come from the MILAGRO experiment, which 
is a wide angle water Cerenkov detector with excellent angular resolution [1 1]. In 
view of the importance of galactic diffuse gamma-ray measurements at TeV ener- 
gies and the availability of several imaging Cerenkov systems, we have conducted 
a detailed simulations- based study of Cerenkov images recorded in the TACTIC 
telescope to identify additional non-imaging parameters which have the potential 
to provide clear segregation between events initiated by isotropic gamma-rays and 
those of cosmic ray hadron origin. Here we report our success in segregating these 
two types of events, recorded by an imaging telescope, by exploiting the fractal 
and wavelet characteristics of the image and use our results to obtain a preliminary 
estimate of the galactic diffuse gamma-ray flux at > 1 TeV energy from a 36h long 
TACTIC data set. 



2 TACTIC system 

The TACTIC ( TeV Atmospheric Cerenkov Telescope with Imaging Camera ) 
gamma-ray telescope is operating at Mt. Abu ( 24.6° N,72.7° E, 1300 m asl), In- 
dia, for the last several years[12] . It essentially comprises a light collector (F/l 
type) of ~ 9.5m 2 area which is made up of 34 x 0.6m diameter,front coated spher- 
ical glass mirrors prealigned to produce an on-axis spot size of 0.3° diameter at 
the focal plane. The telescope is equipped with a 349-pixel photomultiplier tube- 
based imaging camera with a field of view of 6° x 6° and uniform pixel resolution 
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of~ 0.31 . The inner most 225(15x15) pixels are used for generating the event 
trigger based on the 3NCT (Nearest Neighbour, Non-Collinear Triplet) topological 
logic[13] . Operating at a 7-ray threshold energy of ~ 1 TeV, the telescope records 
the Cerenkov images of gamma-ray and cosmic ray hadron-initiated showers with 
a typical rate of ~ 2Hz in the vertical direction. The imaging data is subjected 
to standard image cleaning procedure before being parameterized in terms of the 
standard image shape and orientation parameters. A detailed description of TACTIC 
telescope and the data reduction procedures is available in [12]. 



3 Simulations 



The CORSIKA (version 5.61) air shower simulation code, developed by the KAS- 
CADE group [14- 16] has been used for simulating the extensive air shower develop- 
ment and production of Cerenkov photons in the atmosphere at Mt. Abu by primary 
cosmic-ray protons and gamma-rays. A total of 6000 showers each have been gen- 
erated for cosmic gamma-ray primaries with energies 2, 5, 10, 20, 30 and 50 TeV 
and an equal number for cosmic ray protons with energies 4, 10, 20, 40,60 and 100 
TeV respectively. The primary particles in both the cases have been considered to 
be isotropically incident, within the trigger field of view of 4°x 4° of the TAC- 
TIC telescope, from a direction making an angle of 50° with the vertical. The high 
value of zenith angle has been deliberately chosen to take advantage of the larger 
effective collection area(alongwith a proportionate increase in the threshold energy) 
available in the large zenith angle mode of operation[17]. Atmospheric extinction 
of Cerenkov photons has been taken into account using appropriate wavelength- 
dependent extinction coefficients and the photon distribution in the 349-pixel cam- 
era derived after appropriate ray tracing of the photons for all events passing the 
3NCT trigger criterion[18] . An additional data set of atmospheric Cerenkov events 
initiated by primary gamma-rays of 2 TeV energy, incident from a zenith direction 
of 20°along the principal axis of the TACTIC reflector has been compiled for com- 
parison with TACTIC on -source observations. The derived Cerenkov images have 
been subjected to the usual flat- fielding, image cleaning and calibration procedures 
before being subjected to the fractal and wavelet analysis as described below. 



4 Comparison with TACTIC data 



In order to check whether the simulations reproduce exactly the expected behaviour 
of the TACTIC telescope, we have compared the imaging analysis results with re- 
sults from actual on and off source data scans respectively.The image parameters 
have been calculated as per the procedure described in[19] . It needs to be kept 
in mind that a predominant majority of Cerenkov events recorded by the TACTIC 
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telescope in both the on-and off source scans is of cosmic ray origin(mainly pro- 
tons)while the on-source scan may also contain a possible small gamma-ray signal 
from the source which is being tracked by the telescope. Fig.l shows a compari- 
son of the distribution of image parameters L,W,D and Alpha for simulated 4 TeV 
energy cosmic ray proton initiated events and an equal number of events from a 
TACTIC off-source scan recorded during a Crab nebula observation compaignAs 
expected , the image parameter distributions for the two samples completely over- 
lap each other showing that the simulations are able to reproduce the expected 
behaviour for cosmic ray proton events. This conclusion is reinforced by the results 
presented in Fig. 2 where we show a similar comparison between gamma-ray like 
events selected by using the prescribed imaging cuts to a TACTIC on -source data 
set comprising 259619 events and the results from an equal number of simulated 
2 TeV gamma-ray events incident along the principal axis of the light collector in- 
clined at an angle of 20° with respect to the vertical[20]. Since the majority of the 
events in the experimental data set again result from cosmic ray protons,erroneously 
classified as gamma-ray like events from their image parameter values,we expect- 
edly find a much poorer agreement between the two distributions even though the 
Crab on-source data contains a statistically significant gamma-ray signal as re- 
vealed in the Alpha distribution plot. The above tests conducted with the on and 
off source data sets and the simulated data showing that the simulations are able 
to reproduce exactly the expected behaviour of the TACTIC telescope,thus validate 
the use of the simulated data sets for search of other non-imaging parameters for 
gamma/hadron segregation as described below. 



5 Fractal and Wavelet analysis of Cerenkov images 

The pattern recognition technique for gamma-hadron separation employs fractal 
and wavelet analysis of Cerenkov images to exploit differences in the structure 
of the recorded images due to differences in the longitudinal and lateral develop- 
ment of showers initiated by the two primary species. Showers initiated by pri- 
mary gamma-rays and protons produce different numbers of relativistic electrons 
and muons and also have different heights (and lateral extents) of shower maxi- 
mum. These differences lead to structures in the Cerenkov image on different scale 
lengths which can be parameterized through multifractal and wavelet dimensions 
and moments of different order[21-23] . The usual procedure for calculating these 
parameters involves dividing the Cerenkov image into M = 4, 16, 64, 256 equally 
sized, non-overlapping cells and finding the number of ADC counts in each cell. 
The multifractal moments of order q(q = 2— 6) are calculated from the relation 
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Fig. 1. Comparison of image parameter distribution (L,W,D and Alpha) for simulated 4 TeV 
proton-initiated events and actual TACTIC data recorded during a Crab off-source scan 

where N is the total number of ADC counts in the image, Kj is the ADC count in 
the ) th cell in particular scale and Kj + i in the ) th cell in the consecutive scale and 
q= 2, 3, 4, 5, 6 is the order of fractal moment. The wavelet moments are obtained 
from 
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Wavelet moment W q (M) = ^ 
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The fractal scale length is given by v = log 2 M 

Both the fractal moment G q as well as the wavelet moment W q exhibit a propor- 
tionality with the scale length which can be expressed as, 



G q (M) oc M Tq and W q (M) oc 



(3) 



The fractal and wavelet dimensions of order q are then determined from the rela- 
tions 
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Fractal dimension D q 



q-1 



(4) 



Wavelet dimension B q 
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where r q and (3 q are slopes of the line drawn between natural logarithm of G q (M) 
and W^M) with scale length v respectively. Since wavelet dimensions are more 
sensitive to local structures in the photon distribution in Cerenkov images, they are 
expected to have better efficiency for gamma-hadron separation as proton-initiated 
showers are known to have a non-uniform structure due to contributions from pions 
and muons[17]. 



6 Results from simulations 

Fig. 3 and Fig. 4 show a comparison of the frequency distributions of fractal and 
wavelet dimensions D 2 ,B 2 ,D 6 and B 6 for simulated Cerenkov images produced by 
isotropic( within TACTIC field- of- view) 2 TeV gamma and 4 TeV proton- initiated 
atmospheric Cerenkov events. Also shown for comparison are the B 2 ,D 2 and B 6 ,D 6 
distributions for a TACTIC off- source data set of thirty six hour duration. For q> 
6 both the fractal and the wavelet moments exhibit saturation as a function of scale 
length so that the fractal and wavelet dimensions D 2 , D 6 ,B 2 and B 6 represent the 
two extreme values of these parameters which can be tested for gamma-hadron 
separation.lt has been pointed out earlier[22]that the fractal dimensions D 2 and D e 
distinguish between images containing a single large peak in the centre and images 
with several peaks or a smooth distribution while the wavelet dimensions examine 
differences in neighbouring pixels at different length scales. 

We find from Fig. 3 and Fig. 4 that in the case of fractal dimensions the frequency 
distributions of D 2 and D 6 for gamma-ray and proton-initiated events are very 
broad and overlap over most of the range except for a small peaked structure at 
D 2 = 0.60 and D 6 = 0.50 for gamma-ray events. However, in the case of wavelet 
dimensions B 2 and B 6 , we find that the frequency distributions for gamma and 
proton-initiated events are well separated from each other. The best segregation is 
observed for the parameter B 6 where the frequency distribution for proton-initiated 
events is found to be strongly peaked at B e = 1 .75 while the frequency distribution 
of S 6 for gamma-initiated events produces a slightly broader peak with a prominent 
maximum at B Q = 5.75. For each of the four parameters we have chosen a domain 
of parameter values which passes the maximum number of gamma-ray initiated 
events and rejects the largest number of cosmic ray proton initiated events to quan- 
tify the gamma-hadron segregation capability. The efficiency for gamma-hadron 
segregation has been quantified in terms of the quality factor Q in each case from 
the relation. 
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Fig. 2. Comparison of image parameter distribution (L,W,D and Alpha) for simulated 2 TeV 
gamma-ray initiated events(incident along the reflector axis from 20° zenith direction) with 
actual gamma-ray like events extracted from a TACTIC on-source scan of Crab nebula. 
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where A^ 7 is the number of gamma ray events accepted out of the total number 
(N jT ) of gamma-initiated events after applying the parameter cut. A^represents 
the number of CR proton events retained out of the total (N pT ) proton-initiated 
events after applying the parameter domain cut. The highest quality factor (Q=8) 
is obtained in the case of B 6 , suggesting that this parameter can provide > 99% 
gamma-hadron segregation. As can also be seen visually from Fig. 4, > 90% pro- 
ton events have B % < 2, while > 90% gamma-ray initiated events have B 6 > 5, 
indicating the high efficiency for gamma-hadron segregation using this parameter. 
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Fig. 3. Comparison of distributions for Fractal and Wavelet dimension parameters D2 and 
B2 for simulated 4 TeV proton initiated events and 2 TeV isotropic gamma-ray initiated 
events. The vertical lines represent the parameter cut-off value used to estimate the gam- 
ma-proton segregation efficiency(~ Q). For details refer to the text. 

7 Experimental results 



The TACTIC telescope has been observing the Crab nebula and a few extragalactic 
sources regularly since 2000 in the on-source mode. Occasionally in order to obtain 
a realistic estimate of the background, the telescope is used to track an off- source 
region which is a region of the sky close (~ few degrees away in declination but 
at the same right ascension)to the candidate source and is known to be devoid of 
any TeV gamma-ray source.We have subjected several Crab off-source data sets of 
~ 36h duration in all to the above fractal and wavelet analysis to isolate the back- 
ground diffuse gamma-ray events(if any) from the recorded events using the B 6 
parameter cut. Fig. 5 shows the frequency distribution of B 6 for the events recorded 
during such off-source scans from a zenith direction of ~ 50° with 3h < a < 7h 
and 5 = 22° .OOmin. 52 sec as off-source coordinates. The total number of events 
recorded during the Crab off- source scans was 259600, corresponding to an event 
rate of ~2.0 Hz, which translates into a threshhold detection energy of 4 TeV for 
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Fig. 4. Comparison of distributions for Fractal and Wavelet dimension parameters Dq 
and Bq for simulated 4 TeV proton initiated events and 2 TeV isotropic gamma-ray ini- 
tiated events and actual data recorded during a 36h off-source scan. The lowest panel 
of figure shows an expanded version of the lower part of Bq distribution. The vertical 
lines represent the parameter cut-off value used to estimate the gamma-proton segregation 
efficiency^ Q). For details refer to the text. 

cosmic ray proton-initiated events (~2 TeV for gamma-ray initiated events) as al- 
ready considered in the simulation exercise. While 220320± 469 events out of the 
total events ~ 85% fall within the proton domain of < 2 ,only 154 ± 13 events 
fall in the gamma-ray domain of > 5. Assuming the gamma domain events to re- 
sult from diffuse galactic gamma-rays and taking into account the exposure time 
(36h) TACTIC effective area (7.07 x 10 8 cm 2 ) and the TACTIC field of view(3.8x 
10~ 3 sr),the upper limit of diffuse gamma ray flux (<p UL ) above 2 TeV within 40% 
systematic error [20]is found to be (2.4 ± 0.3)10~ 10 photons cm~ 2 s _1 sr _1 .This 
upper limit is plotted in Fig. 6. 

Although this preliminary result is based on simulation-based fractal parameter cuts 
which are yet to be optimized, the above estimate for galactic diffuse gamma-ray 
flux above 2 TeV compares favourably with the upper limit of 3 x 10~ 8 cm _2 s _1 sr _1 
above 0.5 TeV reported by the Whipple group[9]and <1.96x 10~ 10 cm _2 s _1 sr _1 
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Fig. 5. Wavelet dimension Bq distribution of Cerenkov images recorded by TACTIC tele- 
scope during ~ 36h off-source scan. 

above 10 TeV reported by the Tibet air shower group[24] . Fig. 6 also shows the mea- 
surements and upper limits obtained by other experiments observing the Galactic 
region in the energy range below and beyond 2 TeV. Analysis of simulated data has 
revealed that fractal and wavelet parameters has the capability of better segregation 
of the cosmic-ray and gamma-ray initiated Cerenkov events. In order to test this 
hypothesis, we used fractal and wavelet analysis approach on one spell of 66h Crab 
nebula on-source observation data collected during Dec2007-Jan08 using TACTIC 
telescope and results obtained here showed better matching of experimental and 
simulated non-image parameters(S 6 ,D 6 )[25]. 



8 Conclusions 



It has been shown on the basis of a fractal and wavelet analysis of simulated at- 
mospheric Cerenkov images, recorded in the TACTIC telescope, that it is possible 
to segregate cosmic ray proton-initiated events from diffuse gamma-ray initiated 
events with high efficiency using the wavelet dimension B 6 as the differentiating 
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Fig. 6. TACTIC > 2 TeV upper limits on the diffuse gamma-ray emission within galactic 
range of |6|< -5° and |Z|< 200° 

parameter. Application of data cuts based on this parameter have enabled us to 
provide a preliminary upper limit estimate of the diffuse galactic gamma-ray flux 
above 2 TeV energy from a 36h off-source data set with the TACTIC telescope. 
We hope to improve the simulations based results by optimizing the parameter cuts 
after more extensive simulations taking into account the cosmic ray spectrum and 
TACTIC characteristics. 
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